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Abstract 

We present an experimentally guided, multi-phasic, multi-species ionic gel model to com¬ 
pare and make qualitative predictions on the rheology of mucus of healthy individuals (Wild 
Type) versus those infected with Cystic Fibrosis. The mixture theory consists of the mucus 
(polymer phase) and water (solvent phase) as well as several different ions: H+, Na'*" and Ca^"*". 
The model is linearized to study the hydration of spherically symmetric mucus gels and cal¬ 
ibrated against the experimental data of mucus diffusivities. Near equilibrium, the linearized 
form of the equation describing the radial size of the gel, reduces to the well-known expression 
used in the kinetic theory of swelling hydrogels. Numerical studies reveal that the Donnan 
potential is the dominating mechanism driving the mucus swelling/deswelling transition. How¬ 
ever, the altered swelling kinetics of the Cystic Fibrosis infected mucus is not merely governed 
by the hydroelectric composition of the swelling media, but also due to the altered movement 
of electrolytes as well as due to the defective properties of the mucin polymer network. 
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1 Introduction 

Mucus is a polyelectrolyte biogel that plays a critical role as a protective, exchange and transport 
medium in the digestive, respiratory and reproductive systems of humans and other vertebrates [|Tl|2l. 
Its swelling mechanism are of special interest because of its role in understanding a variety of 
diseases including cystic fibrosis (CF) [HI HI |5l . The conformation of the long-chain, negatively 
charged mucus glycoproteins depends strongly on factors such as pH, ionic strength and ionic bath 
composition [|6||. Mucin is present in secretory vesicles at very high concentrations where they 
are shielded primarily by a combination of divalent ions (e.g., Ca^"*") [|3. Experiments show that 
the mucus gel may swell explosively, up to 600-fold, in times that are the order of a few seconds, 
a process not observed in hydrogel swelling [|8] |9||. Experiments also confirm that this rapid and 
massive expansion of the mucus gel is driven by an exchange of calcium in the vesicle for a 
monovalent ion (e.g., Na~*’) in the extracellular environment [Q. This is because, calcium being 
divalent, must balance two negative charges rather than one. Hence, a divalent Ca^+ ion can act 
as a ‘cross-linker’ between two polymer strands, allowing much tighter condensation than when 
the negative charges of the network are shielded by monovalent ions (Fig.[^. Eurther, experiments 
demonstrate that the exocytosed mucin is recondensed if the calcium concentration of the ionic 
medium is increased sufficiently [flOll . These observations indicate that the amount and the nature 
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of the salt dissolved in the solvent determines the initial and the equilibrium eonfiguration of these 
gels. 

Experimental findings have reported that (unlike hy¬ 
drogel swelling) a simple osmotic pressure difference of 
solute particles across the mucus gel cannot explain the 
massive and the explosive post-exocytotic swelling obser¬ 
vations m. Mucus hydration results from a balance be¬ 
tween osmotic forces, diffusional hindrance (of the ensem¬ 
ble of entangled polymer mesh of the mucus), polyionic 
charges of the mucin and the other fixed polyions en¬ 
trapped within the gel’s matrix, the last two factors being 
the main idea behind the forces via Donnan potential IfT^ . 

The Donnan potential is not constant, but modulated by 
the concentration of free cations and polycations and the 
pH of the hydrating fluid [fT3ll . Marriott et al. observed 
that a decreased monovalent ion concentration, or an in¬ 
creased Ca^"*" concentration, in the airway surface liquid 
was a limiting factor in mucus hydration in CF-infected 
mucus llT4l . Lisle indicated the role of defective process¬ 
ing of mucin polymers, via increased sulfation, in the abnormal mucus hydration in CF flTSl . In¬ 
creased sulfation and sialilation found in CF mucins ffTbl are important for mucus swelling, since 
the high affinity of sulfate residues for calcium drastically shifts the normal monovalent/divalent 
ion-exchange properties of mucins and their swelling characteristics following the exocytotic re¬ 
lease [fTTIl . In summary, besides the ionic composition of water on the surface of the airway, the 
defective processing of mucins and its abnormal ion exchange properties are essential factors to 
explain the characteristically deficient hydration and rheology of mucus in CF. However, the precise 
quantitative link between these elements and mucus hydration still remains unknown [fTFI . 

The theory to understand the swelling and deswelling of ionic gels has a long history beginning 
with the classical work of Flory [fT^ [20l l2Tll and Katchalsky [|22| (see also [|2^ l24l l. An early 
study to understand the kinetics (and not simply the equilibria) of swelling and deswelling, was 
done by Tanaka and colleagues [|^ l26l who developed a kinetic theory of swelling gels viewing 
a gel as a linear elastic solid immersed in a viscous fluid. Although they neglected the motion of 
the fluid solvent, the model reasonably explained the swelling of a gel to its equilibrium volume 
fraction. Their work gave rise to the concept of gel diffusivity as a way to characterize the kinetics 
of swelling. The gel diffusivity is defined as D = where L is the equilibrium size (length 
or radius) of a gel and r is the time constant of exponential swelling toward the equilibrium size. 
They found that expansion of the gel was governed approximately by a diffusion equation with 
diffusion coefficient D. However, how the diffusivity of the gel, or more generally the kinetics of 
swelling (perhaps with large changes in volume fraction) is affected by the movement of the ions 
and defective properties of the mucus polymer network is poorly understood. 

Subsequent studies relaxed the assumption of linear elasticity by defining the force on the gel 
to be the functional derivative of the free energy function for the polymer mesh [|T71 [28l [29l [30l 
[ 3 TI . However, most of these works neglected the fluid flow that must accompany swelling. Wang 
et al. [l32l added fluid flow by application of two-phase flow theory but considered only small 
polymer volume fractions and small gradients in the volume fraction. Burning and Morman llTTll 



Figure 1: Schematics showing how the 
ion-displacement changes between cal¬ 
cium ion (Ca^’*') and a monovalent ion 
(e.g., Na+) leading to changes in the net¬ 
work structure of the mucus matrix and 
its eventual expansion. 
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also used continuity equations to describe the flow of solvent and solution in the gel, but used a 
diffusion approximation with a constant diffusion coefficient to determine the fluid motion. More 
recently, Wolgemuth et al. |[T3l extended these theories to study the swelling of a polyelectrolyte gel. 
However, these current state of the art models do not couple the binding/unbinding of the ions to 
the network micro-structure, an idea which is critical in describing the bulk mechanical properties 
of the polymer [34]| . 

The purpose of this paper is to provide a new comprehensive model detailing the swelling 
kinetics of mucin-like ionic gels and use it to calibrate the rheological data of spherically symmetric 
mucin granules, exocytosed from the goblet cells of CF-patients versus those released from the 
non-CF individuals. The novel feature of this model is the coupling of the dynamical motion of 
the swelling gel with the ionic binding to the polymer network. We use this model to show how 
the altered kinetics of the CF-infected mucus depends on the electro-chemical environment of 
the swelling media as well as the rheological properties of the polymer network. The next section 
presents the details of this model (Q, including equations of motion ( §2.1[ ) and chemistry of binding 
reactions (f2.2). Sections ^2.3 and §2.4| outline the linearized analysis of spherically symmetric 
swelling gels and methods to estimate model parameters, respectively. The results pertaining to 
the swelling kinetics and the equilibrium configuration of these ionic gels under different chemical 
stimuli are presented in ^ We conclude with a brief discussion of the implication of these results 
on the pathophysiology of the infected mucus. 


2 Multi-species, multi-phase mucus-gel model 

The polymer gel is modeled as a multi-component material, consisting of k different types of parti¬ 
cles. Specifically, this material consists of solvent particles, polymers, and several small molecular 
ion species. The polymer is assumed to be made up of monomers (i.e., charge units), denoted 
as M, each of which carries a single negative charge. The positively charged ions in the solvent are 
Hydrogen (H+), Sodium (Na"^) and Calcium (Ca^+). The negatively charged ions could include 
Hydronium (OH“) and Chloride (Cl“). Because the negatively charged ions are assumed to be not 
involved in any binding reactions with the gel, acting only as counterions to positive charges, we 
identify these ions by the name chloride. The binding reactions of the positively charged ions with 
the monomers are 

M- + H+ ^ HM, M- + Na+ ^ NaM, M” + Ca^+ ^ MCa+, M” + MCa+ ^ MaCa, 

k—ri k_c k_x 

( 1 ) 

where k^; and k.^, for C = h,n, c, x, are the binding and the unbinding rates, respectively. We 
assume that all the binding sites/charge sites are identical although the binding affinities for the 
different ions are different. The species M 2 Ca are cross-linked monomer pairs, and the species 
M“, MCa'*", NaM and HM are different monomer species, all of which move with the same poly¬ 
mer velocity. The ion species are freely diffusible, but because they are ions, their movement is 
restricted by the requirement to maintain electroneutrality. Finally, because a small amount of water 
dissociates into hydrogen and hydronium, we are guaranteed that there are always some positive 
and negative ions in the solvent. 
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2.1 Equations of motion and interface conditions 

Suppose we have some volume 1/ of a mixture comprised of k types of particles (including polymer, 
solvent and ion species) each with particle density (number of particles per unit volume), nj (x, t) 
and particle volume, Vj, moving with velocity, Vj(x, f), j = 1,..., k. The subscripts (1, 2) denote 
the polymer and the solvent phase, respectively. We denote these phases with subcripts (p, s) respec¬ 
tively. The volume fractions for polymer and solvent are 9p = VpUp and Og = Vsns- Conservation of 
polymer implies 

riR 

^ + V • (v/p) = 0. (2) 

We assume that the other molecular species do not contribute significantly to the volume. Therefore, 
the volume fractions, 9p + Og = 1. It follows from Q and a similar conservation argument for 
solvent that 

V ■ {9gyg + 9pWp) = 0. (3) 

The motion of the polymer and solvent phase of this multi-component mixture is governed by the 
Stokes equation for Newtonian fluid, which are 

V ■ (6'pcrp(vp)) - 0p(vp - v^) - —Vpp = 0, (4) 

Pg Pp 

and 

Q Q _ ^ Q 

V ■ (^^^^(Vs)) - i—<^p{ys -Vp) -4 - -V/Ug = 0, (5) 

J>3 

where crj(v) = ^(Vv-f Vv^) -f Aj/V • v, are the viscous stresses, rjj > 0 and Xj are the viscosities, 
^ is the drag coefficients, (pp = ‘^plnp+ns, (pg = '^“/up+ns, (pj = for j > 3 (assuming that 

the ions are dissolved in the solvent), are the polymer, solvent and ion species per total solvent 
particle fractions, respectively. The first and second terms in PDE represents the viscous 
forces and drag forces due to the friction between the two phases, respectively. The third term in 
PDE Q represents the force due to the dissolved counter-ions in the solvent (osmotic effect), while 
the last terms in PDE are the chemical forces due to the respective phases. The chemical 
potentials are given by 


/ip = ksTMp + Zm^e + PpP, 

IIg = ksTMg -f PgP, (6) 

where Mp and Mg represent the entropic contribution to the respective chemical potentials (de¬ 
scribed later in this section), ks is the Boltzmann constant, T is temperature, P is pressure, Zm is 
the average charge per monomer (defined later in §2.2[ ), is the electric potential arising from the 
unbalanced charges. For the ion species, the particle volume, pj, is effectively zero so that 

= kBT{\n(pj+ 1-2ai) + Zj^e, j > 3, (7) 

where Zj is the charge on the jth ionic species. In calculating particle fractions, (pj, we assume that 
the particle density of the ions is insignificantly small compared to the particle density of polymer 
and solvent, ^Ug + Up. The ion species satisfy the force balance 

ijUjiwg - Wj) - rijVnj = 0, j > 3. (8) 
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Finally, since there is a free moving-edge to the gel, on one side of whieh (inside the gel) 9~ = 6p, 
and on the other side of whieh (outside the gel) 9^ = 0, = 1, implying that there is no polymer 

outside the gel. The superseripts (-I-, —) denote regions inside and outside the gel, respeetively. The 
interfaee eonditions are 


ksT 




n = 


+ Zm^e + 


ksT 


n. 


(9) 


and 


where 


)) n = ^ -aj + aJ - n, 

0s + ^ (^0s + /^o) 5 

Ms = ln0, + (/.p - a/ + ^ . 


( 10 ) 


( 11 ) 


N is the number of monomers per polymer-ehain and Tq is a referenee temperature. The hydrostatie 
pressure is P~ = P, P^ = 0 and the normalized eleetrostatie potential (T'e = '^‘^/kgr) is, = T'e, 
'1/0 = 0. The normal to the free surfaee is denoted by n. The term aj in the solvent ehemieal 
potential (Eq. (Ill) is the total ion partiele fraetion (a/ = represents osmotie 

pressure as eharaeterized by van’t Hoff’s law. The quantities x, /^o, A^o Flory interaetion 

parameter, standard free energies for pure polymer and pure solvent respeetively. These are related 
to the eross-link fraetion, a (the fraetion of monomers bound with ealeium). 


X — ^(^1 + 62) — 2 — —j ei — eiCK, 

/^o = + 2^3®, 

/^o = -e2|. (12) 

ei {i = 1,..., 4) are the nearest neighbor interaetion energies parameters of the various monomer- 
monomer and monomer-solvent pairs and 2 ; is the number of interaetion sites on a lattiee (i.e., 
eoordination number) [1^ . Eliminating P from Eqns. we find a single interfaee eondition 

+ f^s(v0)) n = Snetn, (13) 


where Snet is the net swelling pressure at the interfaee. 


^net _ ^ ^ I ^ 

7 rri rri ' 6 ' 

ksT Um IZs T Us Um 


(14) 


In the above equation, the fourth term represents swelling pressure eoming from any eleetrie eharge 
on the monomers, referred as the Donnan pressure {Zp^e is the eorresponding Donnan potential) 
whereas the last two terms represents the osmotie swelling pressure eoming from the differenee 
between the eoneentrations of ions dissolved in the gel and those dissolved in the bath. Sirear et al. 
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[|^ derived these equations of motion and the interfaee eonditions using the standard variational 
arguments to minimize the rate of work dissipated within the polymer and the solvent. Unlike 
previous theories, this model aeeurately eaptures how the binding/unbinding of the dissolved ions 
infiuenees the motion of the swelling gel. The ehemistry of the dissolved ions and the eharged 
polymer speeies is deseribed next. 


2.2 Ionization chemistry 

This seetion formulates a mathematieal model to represent how the ehemieal speeies move and 
reaet. Let the eoneentrations per total volume of the polymer speeies be denoted by x = [M 2 Ca], 
m = [M“], V = [NaM], w = [MCa+] and y = [HM], with the total monomer eoneentration 


rriT = m + 2x + w + v + y. 


(15) 


The eoneentrations per solvent volume of the ion speeies are denoted as c = [Ca^’*'], n = [Na"*"], 
h = [H+], and ci = [Cl~]. With eoneentrations expressed in units of moles per liter, the relationship 
between ion partiele fraetions (j)j and eoneentrations Cj is (f>j = UsNaCj, where Na is Avagadro’s 
number. To deseribe the ehemieal reaetions, we use the law of mass aetion. Since all the monomer 
species are advected with the polymer velocity Vp, the monomer species evolve according to 

+ V ■ (vpj) = - RJ, j = X, x, w, y (16) 

where Rj are the forward and backward rates for the monomer binding reactions in Eqn. (fll), 
respectively. The monomer concentration, m, is obtained from Eqn. ( jlsj ) and rriT = ^p/v^Na. Under 
the assumption of fast chemistry, we set the right hand side of the PDE ( [T^ to zero, which reduces 
into the following set of equations for each of the monomer species 


w = 




me, 


V = 




mn, y = 


K. 


mh, 


X = 




m^c, 


(17) 


where = AKc- The above expressions ( [Tt] ), 

assume that the unbinding (dissociation) reactions are ionization reactions that require two “units” 
of solvent. We take the unbinding reaction rates to be for C = c,x,n,h and because calcium 

is a divalent ion, 2kx = kc and k-x = 2k-c. 

Similarly, under the assumption of fast diffusion and chemistry, the law of mass action for the 
ion species reduces to 

C = (18) 


with symbol C = c,h,n,ci, Zn = Zh = = 2 and z^ = —1. The subscript ‘b’ denotes 

the corresponding bath concentrations [|35l. The electrostatic potential, is determined by the 
electroneutrality constraint inside the gel, namely. 


{2c + n + h- ci)6s + ZmmT = 0, 


(19) 


where Zm is the average residual charge of the unbound monomers which depends on the amount 
of binding with ions, 

ZmmT = w — m. (20) 
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Since both the eleetrostatie potential and polymer partiele fraetion are assumed to be zero outside 
the gel, eleetroneutrality in the bath requires that 

2cb + nb + h-cih = 0. ( 21 ) 


Finally, the erosslink fraetion (or the fraetion of monomers bound with ealeium) is a = where 
Eqns. ( [T5|[T7] ) define the eoneentrations, m-r, x, respeetively. 

In summary, the dynamieal motion of a freely swelling mueus-gel is modeled by the system of 
equations ineluding the mass eonservation PDE Q, total volume eonservation PDE Q together 
with foree balanee P D E S (|4]j^ and interfaee eondition Eqn. ( [T3] ), subjeet to the eonstraints Eqn. ( [T5| ) 
(monomer eonservation), Eqn. ([T8]) (ion motion) and Eqn. ([T9|) (eleetroneutrality). 


2.3 Linearized analysis of spherically symmetric swelling gels 


We now eonsider the swelling kineties of a radially symmetrie mueus gel using a linearized analysis 
of the governing equations, deseribed in the previous two seetions. The volume fraetion Op is non¬ 
zero on the domain 0 < r < R{t), with R (maximum radius of the sphere) a funetion of time. Sinee 
only the finite solutions of the governing equations are of interest (in partieular, finite solutions of 
the foree balanee, pdes inside this domain, we assume that Vp = Vg = 0 at r = 0. The 

volume eonservation, PDE d3Mmpfies the eonstraint OpVp + OgVs = 0 throughout the domain. Using 
this eonstraint, we reduee the two foree balanee equations into a single equation (by multiplying 
PDE Q by volume fraetion Og and PDE (|^ by volume fraetion Op and subtraeting) 


0,V- 




Ap c) 2 

H —-—r v 

p2 Qp 




or \ Og J or \ Og ) 


+ + OpOgV (S„et) = 0, 

VpJ 

whereas the interfaee eondition (Eqn. (fTS])) reduees to 

„ a /gpVpA Ap d 2 A, d f r^0pWp \ 

p dr ^ dr \ Og J r'^ dr ^ dr \ Og J 


at r = R(t). 


( 22 ) 


(23) 


Note that V(^ — (^) = V(Snet), where Snet is the interfaee swelling pressure defined in Eqn. (14 1 . 
The veloeities inside the gel are v“^^ = Vp/g and outside the gel are v^^^ = 0. Next, we assume that 
the linearized veloeity Vp is small (the equilibrium veloeity v* = 0) and that Op = 0* + 50p, where 
0* is the equilibrium polymer volume fraetion. Sirear et al. [[35l gives details on the equilibrium 
solution. Sinee this is a moving boundary problem, it is eonvenient to map the domain 0 < r < R(t) 
onto the fixed domain 0<y <1 by making the ehange of variables r = R(T)y and t = r. Eurther, 
we seek spaee-time variable separated solutions of the form Vp = a(r)/i(r/) and 60p = a(r)/ 2 (r/). 


Under these assumptions the reduced foree balance, pde (22) is 


fH , 2(?7e + 2Xe)R ^ 

+ A.)!, * + 


2A, 


+ 




_(Ve + K)y‘^ Ve + Xe\iyj 


^ ^) fi = -o:Rj:uo;)fL 


and the interfaee eondition, Eqn. (231, is 


2A 

(v, + Xe)f[ + —h = o:^uo;)f2 
y 


at y = l. 


(24) 


(25) 
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The supercripts (') and (®) denote the derivative of the functions with respect to the variables y and 9, 
respectively. The net shear and bulk viscosities are given by rje = 9*rip + 9*ris and Ae = 9*Xp + 9*Xs. 
(j)* are the equilibrium particle fractions for the polymer and solvent, respectively. For a finite 


solution of Eqn. (24), it is assumed that /i(0) = 0. Further, we assume that f 2 = fi - These assump- 


(?7e+2Ae)/2 
7?e + Ae 


,/9 = 




^ + ^) + eimiuo; 


tions reduce Eqn. (24) into the homogeneous, spherical Bessel differential equation whose solutions 
have the form /i = y'^Jn{/3y), where 7 = | 

and n = The Bessel function: 

several solutions satisfying the condition /i(0) = 0, we choose the solution with the lowest order 
n = rimitij which determines the equilibrium radius of the swelling gel, Rj = R{t —)■ cxd). 


r;e+Ae 

are of the first kind of order n. Since there are 


Rf= 


, 2 I 


7 e T Ae 

,7e + 2Ae 


(26) 


The variable-separable, linearized form of solution (linearized near equilibrium) for mass conserva¬ 
tion, Eqn. ([^, gives 

^ (27, 


a dr 


f2 Rfy^ 


^ch 


where Tch is a constant (to be determined later). The first and the last part of Eqn. (271 implies that 
a(r) = Because the polymer is conserved (neither created nor destroyed) the velocity of 

the moving boundary must be the same as the gel velocity at the boundary, Vp{y = 1), 


c) /? 

^ = 1 ) = 


(28) 


The solution to the above equation is numerically computed via Matlab ODE solver ode 15, with 
trivial initial conditions. In particular, we note that if i? i?/, then /i(l) ~ constant and in this 
case Eqn. (281 can be solved exactly, i.e., i?(r) = Rf{l — This expression is identical to 

the radial expansion of the linearized form given by Tanaka et. al [|25ll in their kinetic theory of 
swelling hydrogels. Einally, the time constant of gel-swelling towards the equilibrium size (Eqn. 
is 

f2Rfy^ 


Rch — 




i,=i 


(29) 


which defines the diffusivity for spherical gels, D = 


2.4 Parameter Estimation 

Two sets of data are used to calibrate the model for the kinetics of spherically symmetric swelling 
mucus gels. The first experiment, conducted by Verdugo et al., measures the diffusivity data {R?j 
versus Tch values) of the swelling kinetics of exocytotic mucin granules from four (human) CF- 
patients and three healthy individuals at pH = 7.2, T = 37°C and [Ca^+]fe = ImM as well as 
[Ca^+Jft = 2.5mM [[8l|. In the second experiment, performed by Kuver et al., kinetic data was 
collected from cultured gallbladder epithelial cells from wild-type and CF-infected mice at pH = 
7.0, T = 37°C, [Ca^+]b = 4mM and [Na+]b = 140mM jSH. The microscopic composition of the 
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mucus in both of these experiments were reported identical. Fig.|^ and Fig.[^ presents the results 
from these two experiments, respectively. 

Table lists the values of the parameters used in our numerical calculations. The constants in 
the model are the monomer particle volume, Vp, the solvent particle volume, h>s, the coordination 
number of the polymer lattice, 2 ;, and the nearest neighbor interaction energies, e* (Eqn. ([T^), shear 
and bulk viscosity coefficients of solvent (i.e., water), rfs, A^, respectively, at reference temperature. 
The undetermined constants are the binding affinities of the various cations with the gel, Kc, 
(introduced in Eqn. ( [T7| )) and the polymer viscosity coefficients, /ip, Ap, and the drag coefficient, ^ 
(Eqn. 0). 


Table 1: Constant parameters common to all the numerical simulations. The reference temperature 
for the viscosity coefficients and solubility parameters is fixed at Tq = 25°C. 


Constants 

Value 

Units 

Source 

Repeat unit per chain (N) 

266 

- 

m 

Molecular volume of mucus (z/p) 

5 X 10-20 

m3 

Da 

Molecular volume of water (ug) 

2 X 10-23 

m3 

m 

Shear viscosity of water (rjg) 

8.88 X 10-^ 

Pa s 

[|37l 

Bulk viscosity of water (A.,) 

2.47 X 10-3 

Pa s 

m 

Hildebrand solubility (5p) 

1.0928 (a = 0), 1.3258 (a = 1) 

MPa3/2 

[EHl 

Hildebrand solubility for water (f)^) 

48.07 

MPa3/2 

m 


The radially symmetrically swelling mucus gel has a 3-D configuration, which suggests that we 
choose the coordination number, z = Q, mimicking the 3-D structure of the polymer lattice [fT9ll . 
The standard free energies, fc^To/ip and fc^To/i® and the energy interaction parameters, e* (Eqn. (12 1 ) 
are found from the Hildebrand solubility data, 6i [[39]l . The values for the solubility data for materials 
mimicking mucus glycoproteins are given in a study by Mimura [[38l . The fully un-crosslinked 
(no Ca^+ binding) and fully cross!i nk ed (Ca^+ bound) states are denoted by a = 0 and a = 1, 
respectively. The standard free energy is the energy of all the interactions between the molecule 
and its neighbors in a pure state that have to be disrupted to remove the molecule from the pure 
state. The relation between the standard free energies and the solubility parameters (Table [T]) are 

— ksTofil = 

-ksTo^l = (30) 


where Vg = 2 x 10 cm^ is the volume of one molecule of water at reference temperature. 
To = 298K. The negative sign in Eqn. (30) indicates that /je^o/ip, < 0, since they are the 


interaction energies. Using the relations in Eqn. (30), these values are fixed at ei = 4.84, 62 = 
3.74, 63 = —13.70, 64 = 0. The reference temperature of the experiments was fixed at Tq = 25°C. 

The volume of a mucin oligomer/multimer chain is calculated by mathematically modeling a 
chain comprising N freely jointed cylindrical segments of length I and width d, where N = L/lis 
the number of Kuhn segments (or effective rigid cylindrical segments which determines the different 
conformations a chain can have), and L is the end-to-end length of the chain. Eor a gel forming 
mucin (e.g., human MUC 5 AC, used in the experimental data to calibrate our model), I ~ 0.03 /rm, 
L 8 fim and N = 8/0.03 266 O HOI. The radius of gyration, Rg, (defined as the average 
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distance from all Kuhn segments to the eenter of mass of the ehain) and the pervade volume, V, 
(i.e., the approximate spherieal volume of a sphere with radius Rg) is [l^ 


Rg ^ 


f.S X ^0.2 X jyO.588 

2^5 ’ 


1/ 


TT 


Ax - X R't. 
3 ® 


(31) 


Substituting values, the volume V (= Up) 0.05 /rm^. 

The undetermined parameters, namely the binding affinities K^., Kn, Kc and the rheologieal 
eoeffieients rjp, Xp, ^ are eomputed by minimizing a nonlinear least-square differenee funetion be¬ 
tween the experimental values of diffusivity data {Rj versus Tch) and the eorresponding model 
output (Eqns. ( [2^[29l )), implemented via the M ATL AB non-linear least-square minimization fune¬ 
tion Isqnonlin. These values are found as logio(-^n) = —2.27, log]^o(-^/i) = “3.65, log]^o(-^c) = 
—3.12, T]p = 0.11, Ap = 0.31,^ = 1.97 (wild-type mueus) and logio(-^n) = “2.55, log]^o(-^/i) = 
—3.98, logio(-^c) = “7.12, rjp = 1.02, Xp = 2.87, ^ = 0.21 (CF-infeeted mueus). The eloseness of 
fit between the model (highlighted by the solid lines) and the experimental data points is shown 
in Fig The error bars represent the maximum and minimum deviation from the sample points 
and set at 5% margin of error. Notiee the linear relationship between Rj and Tch, predieted by the 
experiments and aeeurately eaptured by the model. 


3 Results and discussion 

The main idea behind this work is to provide an objeetive eomparison of the swelling properties of 
normal versus unhealthy mueus, immersed in an extraeellular medium with ehemieally eontrolled 
eomposition. Well established results exist in literature whieh detail the relationship between the 
swelling rate, final size of the gel and the swelling time [|25l |2^ . However, these results do not 
explieitly outline how these and other physioehemieal elements influeneing the mueus hydration and 
rheology, depend on the nature of the swelling media. Henee, in §3.1[ we explore the relationship 
between the radial size of the swelling mueus gel and the rheology of the mueus polymer, name the 
drag and the viseosity eoeffieients. The effeet of the ealeium and the sodium ions in the solvent on 
the equilibrium size of the mueus blob are detailed in §3.2|and §3.3[ respeetively. 


3.1 Swelling kinetics 


Numerieal simulations were performed to outline the differenees between the swelling kineties for 
WT and CF-infeeted mueus, by altering the eleetro-ehemieal eomposition of the swelling media. 
Fig. [^presents the radius of the swelling gel, R{t), versus time, r, for WT mueus with ealeium 
bath eoneentrations, C;, = 1 mM and Cfe = 2.5 mM (the ‘dash-dot’ and ‘solid’ eurve, respeetively), 
as well as for CF-infeeted mueus (highlighted by the ‘dotted’ and ‘dashed’ eurve, respeetively). 
The sodium ion eoneentration and the pH in the bath are fixed at = 0 and pH = 7.2. These 
eoneentrations eorrespond to the in vitro eonditions of Verdugo’s experiments dU. 

In partieular, note that the expansion of the CF-infeeted mueus (‘dashed’ and ‘dotted’ eurves, 
Fig.[^ is appreeiably slower than WT mueus (‘dash-dot’ and ‘solid’ eurves, Fig.[^. This differenee 
in swelling profiles is a eonsequenee of small drag eoeffieient of the CF-infeeted mueus gels, relative 


to the viseosity eoeffieients (i.e., the ratio 6 = 


f)Rj]p+^p)+f)p(vs+^s) 


is small for CF-infeeted mueus). 


Using the rheologieal parameters, estimated in (2.4 (i.e., Pp, Xp, 0, the solvent parameters listed in 
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Figure 2: Experimental data of the final size squared versus characteristic swelling time {Tch) 
for the two in vitro electro-chemical conditions mentioned in §2.4. For WT mucus the data points 
are represented by (•) and for CF-infected mucus these points are shown by (■). The temperature 
of the experiments was set at 37°C. The model output (i.e., solid lines) closely fits the experimental 
data points. Numerical simulations are done for the data represented by the points (i) and (ii). These 
are shown in Fig. 
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Figure 3: Swelling radius, R{t), of a spherically symmetric mucin gel versus time, r, with kinetics 
governed by Eqn. (28), for WT-mucus at calcium concentration Cb = 1.0 mM (dash-dot curve), 
CF-mucus at C;, = 1.0 mM (dotted curve), WT-mucus at Cb = 2.5 mM (solid curve), CF-mucus at 
Cb = 2.5 mM (dashed curve). The {Rf, Tch) values for the dash-dot and the dotted curve correspond 
to the data represented by points (i) and (ii) in Fig. [^, respectively. 


Table [^(i.e., As) and the equilibrium polymer volume fraction, 6* (i.e., by solving Eqns. ( [T3| 
191, values shown in Eig. |^); we find that 5 = 82.8655, 78.7264 for WT mucus for concentrations 
Cb = 1 mM and Cb = 2.5 mM, respectively while the corresponding values are 5 = 0.0885, 0.0914 
for CE-infected mucus, respectively. Gels with higher viscosities expand at a slower rate and this 
explains a relatively slow or insufficient hydration due to the defective rheology of the CE-infected 
mucus. Previous numerical studies by Sircar at al. corroborate these results iHTl . 

The corresponding diffusivity data for the curves in Eig. are given in Table The diffusivity 
decreases by 33% for normal mucus and 88% for CE-infected mucus, when the Ca^+ concentration 
in the swelling medium is changed from 1 to 2.5 mM, respectively. Eurther, there is a 76% reduction 
in the diffusivity of CE-infected mucus versus the WT mucus in a solvent containing 1 mM Ca^+, 
in concurrence with Verdugo’s swelling experiments ||8l. A lower diffusivity of the defective mucus 
is a consequence of the slow rate of expansion of these gels (or a bigger time constant, Tch), which 
again lead to the conclusion that the defective rheology of mucus plays a crucial role in the eventual 
mucus hydration. 


Table 2: Diffusivity data for WT and CE-infected mucus at calcium bath concentrations C;, = 1 mM 
and Cf, = 2.5 mM. Rf,Tch and D are measured in /rm, sec and cm^/s, respectively. The sodium 
concentration and the pH in the bath are fixed at N;, = 0 and pH = 7.2. 



1 mM Cd?+ 

2.5 mM Ca2+ 

WT 

CE 

Rf = 2.41, Tch = 0.2, D = 2.90 x 10“^ 
Rf = 2.37, Tch = 0.79, D = 7.09 x 10"® 

Rf = 2.37, Tch = 0.29, D = 1.93 x 10“^ 
Rf = 2.37, Tch = 6.58, D = 0.85 x 10"® 
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3.2 Equilibrium configuration versus bath [Ca^+] 

In the next two seetions, we explore the role of the eleetrolytie eomposition of the swelling media 
on the equilibrium eonfiguration of the healthy versus the diseased mueus. Fig. |^,b,e highlight 
the equilibrium volume fraetion, 6* (whieh determines the swelled/deswelled state of the gel), the 
Donnan potential and the erosslink fraetion, at pH= 7.2 and pH= 5.0, versus the bath eoneentration 
of ealeium, respeetively. 

The emergenee of three different swelling regions are observed in Fig.|^ At low bath eoneentra- 
tions of ealeium (C/, < the gels are immersed in a solvent whieh laeks eations to bind with 

the eharged monomers. This leads to a highly ionized state (i.e., a high average negative eharge per 
monomer) and eonsequently a higher Donnan potential (Fig.|^). Sirear et al. has shown earlier that 
Donnan potential is the dominant meehanism driving the swelling of ionie gels [l35l . Solvents with 
lower pH (swelling profiles represented by the ‘dash-dot’ and ‘dotted’ eurves, Fig.|^) ean furnish 
more H+ ions, leading to a lower Donnan potential whieh translates to a relatively de-swelled state 
(i.e., eomparing the values at pH= 5 versus those at pH= 7.2). 

The hydration properties of the diseased versus healthy mueus are markedly different for small 
to intermediate ealeium eoneentrations, 10“^^ M<Cfe < 10 M (Fig. |^. In partieular, the diseased 
mueus (swelling profiles shown by the ‘dashed’ and ‘dotted’ eurve) shows massive deswelling. The 
explanation for these swelling features in this region is that there is a eomplex interplay of ionization 
via Donnan potential (whieh favors swelling) and the energy gain due to inereased erosslinking 
(whieh promotes deswelling) lf35l . Defeetive mueus gels have higher ealeium binding affinity (i.e., 
log^o-^c = —7.12 for CF-infeeted mueus versus log^g iTc = —3.12 for WT mueus, refer §2.4| 
where these binding affinity values are mentioned). This leads to a relatively de-ionized state in the 
defeetive gels (or a lower Donnan potential, Fig.|^). Consequently, a de-swelled state (mediated 
by the de-ionization due to ealeium erosslinks) dominates. 

At suffieiently high ealeium ion eoneentrations (C?, >10 M) the average eharge on the monomer 
is positive rather than negative. The positive eharge on the monomer is beeause in the divalent ion 
ease, the binding of a monomer with an ion eonverts it from a negatively eharged ion, M~, into a 
positively eharged ion, MCa"^ (Eqn. Q). This has little effeet on the overall swelling pressure sinee 
at high ion eoneentrations the eharge on the gel is small relative to the overall number of available 
ions, and the Donnan potential, vanishes. As a eonsequenee, at suffieiently high C/, the gel (both 
normal as well as the diseased type) behaves as if it is uneharged. In summary, an inereased affinity 
of the CF-infeeted mueus to bind with ealeium (in physiologieally signifieant eoneentration range) 
leads to a highly de-swelled state of the gel. This defeetive binding property of the diseased mueus 
partly explains how the altered movement of the eleetrolytes leads to an insuffieient hydration of 
the gel. 

3.3 Equilibrium configuration versus bath [Na+] 

This seetion reports the salient features of the equilibrium state of the mueus gel immersed in a bath 
eontaining monovalent eations (i.e., Na+, Fig.[^. Again, three different swelling transition regimes 
are seen from these figures: a de-swelled state (i.e., a relatively high polymer volume fraetion, 6*p 
is favored at low (N^ < 10“^ M) and high (N^ > 1M) bath eoneentrations. A higher polymer 
volume fraetion is beeause of a low Donnan potential in these eoneentration ranges, whieh is due 
to either the furnished H+ ions (at low sodium eoneentrations) or exeessive Na+ ions (otherwise). 
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(b) Donnan swelling pressure, Zm^'e 








In small to intermediate sodium eoneentrations (10“^M<Nb < 10”^ M, i.e., approximately in 
the region N;, there is an additional feature with inereasing sodium eoneentration, namely, 

swelling either gradually {pH = 5.0 eurves, Fig.|^) or via phase-transition {pH = 7.2 eurves, Fig. 
[^) followed by de-swelling. The explanation for this swelling feature is a eomplieated interplay 
between ionization via hydrogen unbinding (whieh promotes swelling) and a de-ionization via 
sodium binding (whieh promotes de-swelling). 

However, the volume transition eurves are qualitatively different for the WT and the CF-infeeted 
mueus, e.g., notiee the hysteretie swelling transition for the WT-mueus (solid eurve) versus a double- 
hysteretie eurve for the defeetive mueus (dashed eurve) at pH = 7.2. These differenees, again, arise 
due to different binding affinities (logj^o = —2.27, log^o = —3.65 for CF-infeeted mueus 
versus log^giTn = —2.55, logj^g= —3.98 for WT mueus, values listed in §2.4[ ). Thus, the 
differenees in the ion-binding property of the mueus influenees the not only the swelling-deswelling 
volume transitions but also the equilibrium state of the gel. 


4 Conclusions 

This paper develops a new, eomprehensive, multi-phase, multi-speeies model to quantify the swelling 
/ deswelling meehanism for muein gels. This model explains how the final eonfiguration of these 
gels depends on eomplex interaetions between eompeting effeets that alters the gel ionization, that 
is the Donnan potential, ehanges in the bath eoneentration of ions and their eorresponding binding 
affinity with mueus. Near equilibrium, the radial size of the swelling mueus gel reduees to the 
well known expression for hydrogel swelling [I31 . The diffusivity, D (whieh is similar to the one 
defined by Tanaka’s hydrogel theory) aeeurately eharaeterizes the swelling properties of the muein 
network, ineluding the effeet of ealeium binding on the equilibrium eonfiguration of the gel. The 
reasons behind the experimental outeomes of defieient mueus hydration for CF-infeeted patients 
are examined. In partieular, the defeetive swelling properties of the diseased mueus is the eombined 
effeet of an abnormal mueus rheology (i.e., an abnormally low drag relative to the polymer viseos- 
ity), an altered eleetrolytie eomposition of the extraeellular medium (i.e., an inereased ealeium ion 
eoneentration in the airway surfaee liquid surrounding the mueus) as well as an altered movement 
of ions in the medium (i.e., a higher affinity of the CF-infeeted mueus to bind with ealeium ion). All 
these faetors favor a eondensed state of the gel, resulting in a mueus gel with defeetive mueoeiliary 
transport eontributing to the ehronie airway infeetion found in CF patients. 
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and the temperature are fixed at C?, = 0 and T = 37°C, respectively. 
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